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Entanglement has been recognized as a key resource for quantum computation^ and quan- 



v/-^ I turn cryptography^. For quantum metrology, use of entangled states has been mainly dis- 



cussed as a means of improving the signal-to-noise ratio"* ^ and first demonstration experi 



O 

^ ! ments were performed^. In addition, entangled states were used in experiments for efficient 

p quantum state detection^ and for the measurement of scattering lengths^. In quantum infor- 

+-> ' 

^ . mation processing, manipulation of individual quantum bits allows for the tailored design of 
^3 

R^' specific states that are insensitive against detrimental influences of an environment'. Such 
> . 

^ ■ decoherence-free subspaces^" protect quantum information and yield significantly enhanced 
coherence times^^ Here, we demonstrate precision spectroscopy using a decoherence-free 
subspace with specifically designed entangled states^^ to obtain the electric quadrupole mo- 
ment of Ca+ to be used for frequency standard applications. We find that entangled states 
are not only useful for enhancing the signal-to-noise ratio in frequency measurements but 
that moreover a suitably designed pair of atoms will allow clock measurements in the pres- 
ence of strong technical noise. The applied technique makes explicit use of non-locality as an 



entanglement property and constitutes thus a new paradigm for designed quantum metrol 



Metrology seeks to relate all measurements to a precision determination of frequencies. The 
most precise measurements today are obtained from measuring atomic transition frequencies by 
monitoring the phase evolution of a superposition of the pertaining states. This is achieved by 
applying Ramsey's interferometric technique that can be generalized^ to (maximally) entangled 
states: Under free precession, the two-atom state ^ — ■^{\ui)\u2) + \vi)\v2)) evolves into the state 
^(r) = -^{Wi)W2) + exp(iA^T)|i'i)|i'2)) where the phase evolution rate A<^ = [(-E'ui + ^^2) ~ 
{Ey^ + Ey^)]/ h is proportional to the differences in atomic energies E^^ , Ey^ of the involved levels. 
The real part of the phase factor e*^"*^ is measured by projecting the atoms onto the states \±)k — 

-^{\uk) ±\vk)) and measuring the parity operator P — P{++) +P{ ) — — P{-+) > where -P(±,±) 

denotes the projector onto the state \±)i®\±)2. If^ belongs to a decoherence-free subspace, free 
precession times r of several seconds^^'^^ allow for highly accurate phase estimation. 

In atomic optical frequency standards based on single hydrogen-like trapped ions'^, the tran- 
sition frequency from the S ground state to a metastable Dj state is measured by an optical fre- 
quency comb^^. The Dj state's atomic electric quadrupole moment interacting with residual elec- 
tric quadrupole fields^^ gives rise to frequency shifts of a few Hertz. The shift of the Zeeman 
sublevel \Dj, rrij) in a quadrupole field y, z) — A{x^ + — 2^;^) is given by 



where dE^/dz = iA is the electric field gradient along the potential's symmetry axis z, (3 denotes 

2 



ogy. 



i(j + 1) - 3m/ 
J(2j - 1) 



(3cos2/3- 1), 



(1) 



the angle between the quantization axis z and Q{D,j) expresses the strength of the quadrupole 
moment in terms of a reduced matrix element^^. 

Recently, quadrupole moments have been measured for ^^Sr+, ^^^Hg^ and ^''^Yb+ with a 
precision ranging from about 4% to 12% (refs. ^^"^^). In these single-ion experiments, narrow- 
linewidth lasers are employed to detect electric quadrupole shifts by measuring the transition fre- 
quency from the electronic ground state to the metastable state. Alternatively, the quadrupole shift 
could also be measured by performing Ramsey experiments with an ion prepared in a superposi- 
tion of Zeeman Dj sublevels. Both measurement schemes are subject to phase decoherence. Laser 
frequency noise limits the coherence time if the S state is used as the reference state. Magnetic 
field noise giving rise to first-order Zeeman shifts renders the measurement scheme based on D 
state superpositions difficult. However, both sources of phase noise can be eliminated by replac- 
ing the single ion by a two-ion entangled state prepared in a decoherence-free subspace. Ramsey 
experiments with the two-ion Bell state \l/ = |"^2) + l^^a) |"^4)), where the magnetic quan- 

tum numbers mj of the states |mj) = \Dj, rrii) satisfy nii + m2 = + m^, are not affected by 
fluctuations of the magnetic fields to first order since both parts of the superposition are Zeeman- 
shifted by the same amount (see Fig. [I}?). Also, frequency noise of the laser used for preparing \& 
is relevant only during the comparatively short state preparation and read-out steps but not during 
the long interrogation period. For suitably chosen values of mi, . . . m4, a decoherence-immune 
state can be designed that is sensitive to the electric quadrupole shift (Fig.[T]:). 

The isotope ""^Ca^ has a single valence electron and no hyperfine structure (Fig.[T^). For its 
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(a) (b) m= -5/2 -3/2 -1/2 +1/2 +3/2 +5/2 




Figure 1: Relevant atomic levels of the ^"Ca"*". (a) Bell states are prepared by coherently exciting 
the 5*1/2 ^ -D5/2 quadrupole transition. The 5*1/2 ^ P1/2 is used for quantum state read-out, 
lasers at 866 and 854 serve as repumping and quenching lasers, (b) Energy level shift of the 
|m) = \D^/2,m) states in a weak magnetic field, (c) Energy level shift caused by an electric 
quadrupole potential. Energy levels occurring in the state ^1 = ;^(|o)|«) + | V)|T)) (eq. (O) are 
denoted by the symbols o, v, •, the open (filled) symbols corresponding to levels occupied by 
atom 1 (2), respectively. The combined Zeeman energy of the |o)|«) -states equals the energy of 
the I V) I T) -states, making \E'i immune against magnetic field noise. Since an electric quadrupole 
field shifts the |o)|«)- and | V) |t) -states in opposite directions, \E'i is suitable for a measurement of 
the quadrupole shift. 
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metastable 3^^1)5/2 state (lifetime td^^^ — 1.168(7) s), a quadrupole moment of 1.917eao was 
calculated^^'^°. For the experiment, two "^^Ca^ ions are trapped along the axis of a linear ion trap. 
By applying a static (dc) voltage ranging from 500 to 2000 V to the tip electrodes^"*, axial center- 
of-mass (COM) mode trap frequencies Uz ranging from 850 to 1700 kHz are achieved, leading to 
distances between the ions from 6.2 jjim to 3.9 jjim. A CCD camera images the ions' fluorescence. 
The degeneracy of the Zeeman states is lifted by applying a magnetic bias field of 2.9 G. The ions 
are cooled to the ground state of the axial COM mode by Doppler and sideband cooling^^. For 
coherent quantum state manipulation, the ions are individually excited on the \S) = Si/2{m — 
—1/2) <-> \D) = £)5/2(m = —1/2) transition by a tightly focussed laser beam (laser linewidth 
pa 150 Hz) resonant with either the carrier transition or the upper motional sideband. For a detailed 
description of the experimental setup, see ref.^. 

The Bell state = {\S)\D) ± \D)\S)) / V2 is created with a fidelity of about 90% by 
a sequence of three laser pulses Additional carrier tt pulses transfer the entanglement into the 
D5f2 Zeeman state manifold, thus generating the state 

*i = i=(| - 5/2)1 + 3/2) + I - 1/2)1 - 1/2)). (2) 

This magnetic field-insensitive state is used for a measurement of the quadrupole shift. In the 
vicinity of the trap center, the dc voltage applied to the tip electrodes creates a rotationally sym- 
metric electric quadrupole potential which shifts the energy of the constituents of the superposition 
state ^'i by hAi = 24/5 h6 with respect to each other. Here, hS is the shift that a single ion in state 
I — 5/2) would experience. The two-ion shift is bigger by 12/5 because two ions contribute that 
are prepared in substates which shift in opposite directions. In addition, the presence of a second 
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ion doubles the electric field gradient at the location of the other ion'^. For the shift measurement, 
^1 is let to evolve into ^i(r) = ^(| - 5/2)| + 3/2) + exp(zAir)| - 1/2)| - 1/2)) andcos(Air) 
is measured for r ranging from to 300 ms. Figure shows the resulting parity oscillations^^ at 
a frequency Ai = (27r) 33.35(3) Hz. Disentanglement of ^i(t) by spontaneous decay manifests 
itself as damping of the oscillations with a damping time constant = 587(70) ms close to the 
expected value ^to^^^ = 584 ms. Parity oscillations of the state ^2(t) = -^(13/2)1 — 5/2) + 
exp(iA2r) I — 1/2) | — 1 /2)) are shown in Fig.|2b- Its oscillation frequency A2 = (27r) 36.52(4) Hz 
slightly differs from Ai because of a magnetic field gradient in the direction of the ion crystal 
that gives rise to an additional contribution to the parity signal^^. Both phase oscillation contribu- 
tions can be separately determined by taking the average A = (Ai + A2)/2 and the difference 
Ab' = |Ai — A2I/2 of the signals. By recording parity oscillations at different tip voltages Uups, 
we measure the quadrupole shift A/(27r) as function of the electric field gradient and observe a 
linear change in A{Uups). 

For a precise determination of the quadrupole moment, the dependence of the quadrupole 
shift on the orientation of the magnetic field should be minimized. We investigate the angle depen- 
dence of A by varying the orientation of the magnetic field in a plane that also contains the z-axis 
of the trap. Initialization of the ions in a pure state prior to coherent manipulation is achieved by 
optical pumping on the quadrupole transition (see Methods). Generally, it cannot be excluded that 
the quadrupole field has imperfect rotational symmetry. In this case, the factor 3 cos^(/3) — 1 in 
eq.dU) needs to be replaced by (3 cos^ P — I) — e sin^ P cos(2a) where e characterizes the asym- 
metry and a its direction'^. Note, however, that the additional term vanishes for P = 0. Figure 

6 




Time (s) 

Figure 2: Parity oscillations of the entangled states \E'i and \E'2 at Uups = 750 V tip voltage. 
The experiment is repeated 100 times for each data point. An exponentially damped sinusoidal 
function is fitted to the data, (a) For \E'i, the oscillation frequency is Ai = {2n) 33.35(3) Hz and 
the damping time 587(70) ms. (b) For ^2, the fit yields A2 = {2tt) 36.52(4) Hz and a decay 
time 525(60) ms. For each plot, five different datasets were merged. No data were taken covering 
waiting times around 170 ms. To exclude slow drifts of the initial parity over time, data for waiting 
times up to 20 ms were repeatedly taken. 
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Figure 3: Angular dependence of the quadrupole shift, (a) Quadrupole shift A/(27r) at Uups = 
1000 V as a function of the magnetic field orientation (3 — (3^. The data are fitted by A = + 
Afe cos^(/? — (3q), yielding l3o = 26.9°. The arrow indicates the angle corresponding to /3 = 0. For 
the determination of the quadrupole moment, shift measurements are performed at /3 = since 
errors in [3 enter the calculation only in second order. Note that the error bars are smaller than 
the point size. The inset in the upper right comer shows an enlarged data point, (b) Geometry of 
magnetic field coils. The angle [3 is varied by appropriately changing the current in pairs of coils 
producing the magnetic fields Bhi, Bh2- Field is used to null residual fields in the perpendicular 
direction. 
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[3] shows the sinusoidal variation of A measured at Uups = 1000 V as a function of the angle. At 
the angle that maximizes A, another maximization is performed in the perpendicular direction to 
determine the direction corresponding to /5 = 0. Figure displays the quadrupole shift measured 
at = as a function of the electric field gradient dE^/dz. Using dE^/dz = —mcul/e, the gradi- 
ent is calibrated by a measurement of the ion's oscillation frequency Uz- By fitting a straight line 
A/(27r) = Ao/(27r) + to the data, we determine the slope a = 2.975(2) HzmmW and the 
offset Ao/(27r) = —2.4(1) Hz. The major part of the phase oscillation frequency Aq at ^ = 
is caused by the second-order Zeeman effect that contributes AB2/{2n) = —2.9 Hz at the bias 
field of 2.9 G. The remaining part Aq — is attributed to a residual quadrupole field caused 
by stray charges, ac-Stark shifts are insignificant. The quadrupole moment Q{3d,5/2) = j^ha 
is proportional to the slope a. While the uncertainty in the value of a is smaller than 10^^, the 
main uncertainty in 6(3(i, 5/2) is due to errors in the determination of the angle /3. Assuming 
that the angle can be determined with an accuracy of 3 degrees (see Methods), we determine the 
quadrupole moment to be Q{3d, 5/2) = 1.83(l)eao. 

These results show the viability of the approach of designing specific entangled states for 
precision measurements. For frequency standard applications, minimization of the quadrupole 
shift is imperative^^'^^. Note that, by using designed entangled states, the quadrupole shift could 
be cancelled by averaging the results over states with different magnetic quantum numbers^^. The 
design approach demonstrated here exhibits a new paradigm for metrology: while previously the 
use of entangled states was to enhance the signal-to-noise ratio, the use of a designed pair of 
entangled ions allows an atomic clock measurement to be performed in a subspace free of magnetic 
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Figure 4: Level shifts as a function of the applied external electric field gradient dE^/dz. (a) 
Quadrupole shift. The offset at dE^/dz = is mainly due to the second-order Zeeman effect. The 
quadrupole moment is proportional to the slope Miich is measured with an uncertainty of less than 
0.1%. The inset shows the deviation of the data points from the linear fit. (b) Phase oscillation rate 
AB//(27r) due to the magnetic field gradient calculated from IA2 — Ai|/2. Since the distance d 



field noise. This is achieved with unprecedented precision on an even isotope which opens the 
perspective for optical frequency standards based on a whole variety of easily accessible atomic 
systems not considered so far. 

Methods 

Frequency-resolved optical pumping. Due to the magnetic bias field, the ions can be prepared in 
the \ Si/2-,rn = —1/2) state for an arbitrary orientation of the magnetic field by frequency-resolved 
optical pumping on the quadrupole transition. For this, the \Si/2,m = +1/2) \Dr,/2,Tn = 
—3/2) transition is excited while the repumping lasers on the P3/2 and -D3/2 ^ P1/2 

transitions are switched on at the same time. After an excitation of 1 ms duration, at least 99.9% 
of the population has been pumped to the \Si/2,m = —1/2) state. 

Orientation of the magnetic field. The magnetic field B is set by passing currents Ih^, Ih2, Iv 
through three mutually orthogonal pairs of coils (see Fig. [3]b). We calibrate the field coils and null 
the offset fields by measuring the Zeeman shift of the \Sii2,rn = —1/2) ^ |D5/2,m = —1/2) 
transition as a function of the currents. We estimate that the offset fields can be zeroed to a value 
below 1% of the applied bias field. After nulling the field in the v-direction, we map out the values 
of (//ij, that yield a constant Zeeman shift and fit an ellipse to the data. These current values 
rotate B in a horizontal plane that also contains the symmetry axis (z) of the ion trap. A similar 
procedure is applied for calibrating the v-coil. 

Alignment of the magnetic field with the trap axis. We align the direction of the magnetic field 
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with the trap axis by maximizing the quadrupole shift as a function of the field direction. A residual 
quadrupole potential $s caused by stray charges and by voltages applied to the compensation elec- 
trodes will in general have principal axes that do not coincide with those of the trap's quadrupole 

— * 

potential. This will lead to errors in the alignment of the magnetic field B with the principal axis 
of the trap's electric quadrupole potential ^ups- This is due to the fact that B is aligned with the 
direction Hts that maximizes the quadrupole shift of ^upsi'fT) + ^s(^) instead of the direction Ht 
that would be chosen without $5. Fortunately, $s is much smaller than %ips in the experiment. 
We obtain information about the magnitude of $s by (1) extrapolating the ion oscillation frequency 
LOz{Utips) for Utips — > 0, (2) measuring the electric stray field Es and assuming the magnitude of the 
field gradient to be oc \Es\/r where r is the smallest distance of the ions from the trap electrodes. 
Most importantly, (3) by measuring the residual quadrupole shifts for /3 = and another direction 
with j3 = 26.9°, constraints regarding the shape, orientation and magnitude of $s can be found 
that allow us to calculate the average angle A/3 between fits and fit. We find A/3 = 3° and use this 
value for estimating the uncertainty in the determination of ©(3d, 5/2). 
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